Electrodynamics of a Clean Vortex Lattice 
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We report on a microscopic evaluation of electrodynamic response for the vortex lattice state of a 
model s-wave superconductor. Our calculation accounts self-consistently for both quasiparticle and 
order parameter response and establishes the collective nature of linear response in the clean limit. 
We discuss the effects of homogeneous and inhomogeneous pinning on the optical conductivity and 
the penetration depth, and comment on the relationship between macroscopic and local penetration 
depths. We find unexpected relationships between pinning arrangements and conductivity due to 
the strongly non-local response. 
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The study of vortex electrodynamics in type II super- 
conductors has been ongoing for over 30 years Q. Due 
to the complexity of the problem, however, the theory 
of vortex electrodynamics remains incomplete. For a 
system with perfect translational invariance, it can be 
shown rigorously that the long-wavelength response con- 
sists only of a resonance at the cyclotron frequency 0, 
so that there is no signature of superconductivity in the 
a.c. conductivity. In real systems the response is com- 
plicated by broken translational symmetries associated 
with disorder, pinning defects and even the underlying 
atomic lattice. Relatively few rigorous results exist for 
the general case. One regime which has been exten- 
sively studied is the the dirty limit of vortex core sizes 
£ much larger than quasiparticle mean-free-paths I. In 
this case, the response is approximately local and the 
Bardeen-Stephen ||[| isolated vortex model provides a 
useful picture. Even in the dirty limit, vortex-vortex in- 
teractions play an important role in the electrodynamics 
except at fields very close to H cl and microscopic models 
rapidly become intractable. Instead, one typically in- 
terprets experiments at sub-THz frequencies in terms of 
simpler phenomenological hydrodynamic models ^,^,|llj , 
in which the mixed state is characterized by vortex pin- 
ning, viscosity, and interaction parameters. 

No comparable phenomenology exists for the clean- 
limit vortex lattice, for which I is larger than £ or even 
larger than the distance between vortices. In this case 
the details of disorder within a particular vortex core 
become important and nonlocal response can make 
the connection between single-vortex models || and the 
dense vortex systems tenuous. Interest in the clean limit 
has grown in recent years with the advent of clean sin- 
gle crystals of type II low (s-wave) and high (d-wave) T c 
materials. In particular, the extent to which in-field opti- 
cal conductivity measurements can be used as a probe of 
fundamental properties of high T c materials is a subject 
of great interest @Jl3)- Existing experimental work is 
confusing and possibly contradictory. Resonances have 
been observed in thin films of YBa2Cu307_a (YBCO) 
and were interpreted as being due to vortex core tran- 
sitions [tL2|, the clean-limit manifestation of vortex mo- 



tion [Q. Other, possibly related, resonances observed in 
Bi 2 Sr 2 CaCu208 (BSCCO) [jl3| were interpreted as being 
due to collective vortex motion. In contrast, recent ex- 
periments H in thin films of BSCCO were convincingly 
interpreted in terms of a continuum of extended quasi- 
particle states, ignoring vortex motion. Finally, infrared 
measurements of YBCO crystals have failed to find any 
appreciable field-dependence to the conductivity [jl4| at 
all! 

Much of the theoretical work on electrodynamics in the 
clean limit has been based on one of two extreme pictures; 
in the first quasiparticle states of the vortex lattice 
respond directly to the perturbing electromagnetic field 
without any order parameter response, while in the sec- 
ond (?],[l(| it is the vortices that move, usually rigidly, and 
indirectly disturb the quasiparticle equilibrium. Another 
stimulating approach uses a phenomenological composite 
model to interpret the multiple resonances and magne- 
tooptical activity seen in experiment |l7j . The difficulty 
which arises in any serious calculation, however, is that 
quasiparticles and vortices are not separate entities. A 
time-dependent order parameter, for example one arising 
from vortex motion, creates a time-dependent potential 
for the quasiparticle states from which the order param- 
eter is in turn constructed. 

In this Letter we report on a microscopic calculation 
which treats the order parameter and the quasiparticles 
on equal footing Jl8| , and which naturally captures the 
nonlocal nature of clean-limit vortex electrodynamics. In 
this case, the appropriate language is not that of quasi- 
particles interacting with moving vortices, but is instead 
that of collective modes. The approach is microscopic 
and proceeds in three steps. First we calculate a self- 
consistent solution for the equilibrium vortex lattice us- 
ing the Bogoliubov-deGennes (BdG) equations. Second 
we evaluate the linear response of normal and anomolous 
blocks of the one-particle density matrix to a perturbing 
electromagnetic (EM) field, accounting for vortex motion 
and distortion self-consistently. Third, we use the den- 
sity matrix to find the induced current, and therefore 
the conductivity. Since this is a linear response calcula- 
tion it addresses external fields which produce a vortex 
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displacement that is small compared to the coherence 
length. Because the equations are solved in real space, it 
is a trivial matter to incorporate disorder or other forms 
of broken translational invariance. 

There are two important energy scales in our discus- 
sion. The first is the cyclotron energy hu) c , which is the 
resonant frequency of a translationally invariant super- 
conducting system [||. The second is the quasiparticle 
absorption energy huj qp ~ A 2 /W ~ Huj c (H C 2/H) [ p^[ , 
where A is the s-wave gap, W is the metal band-width, 
and H is the field-strength. This is the predicted en- 
ergy of allowed transitions between vortex core states in 
non-self-consistent calculations which ignore vortex mo- 
tion . Note that the cyclotron frequency drops below 
the quasiparticle bound state transition energy as the 
external field drops below H C 2- Before turning to a dis- 
cussion of our results, we briefly outline the formulation 
of vortex-lattice equilibrium state and response-function 
theory on which they are based. 

We have obtained self-consistent solutions of mean- 
field equations for equilibrium and linearly perturbed 
states of two-dimensional generalized Hubbard model 
Hamiltonians. We use a Nambu vector notation, defin- 
ing C = (ci, c 2 , ... , c 2N ) where, c, = c iT , c. i+N = ctp c\ a 
and Ci a are creation and annihilation operators for site 
i and spin a, and N is the number of lattice sites. The 
mean field hamiltonian can then be written as TL MF = 
c t H BG C where h bg ig the 2N x 2N Bogoliubov-de 

Gennes matrix: 
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Here e, and A are N x N matrices and p, is the chem- 
ical potential. The one-body matrix elements, ej-j can 
be used to model intersite hopping, and diagonal or off- 
diagonal disorder. The magnetic field, which we take 
to be uniform, is accounted for by multiplying all one- 
particle matrix elements by Landau gauge phase factors; 
iji = 6ji exp(iaYi(X j — Xi)) where (JQ, Yi) is the location 
of the i-th lattice site and a = eB/hc. Ay = Vij {cj j.Cjf ) 
where Vij is the Hubbard model interaction strength be- 
tween sites i and j and the angle brackets denote a 
consistently determined thermal average. We enforce 
periodicity |2^,|l| in a supercell containing an even- 
integer number of vortices, typically two. The mean- 
field Hamiltonian can be written in the diagonal form, 
T~t AIF = J2 n Enlnln where E n is an eigenvalue of H BG , 
7„ = J2i Un is a fermionic quasiparticle annihilation 
operator, and U specifies the unitary transformation that 
diagonalizes H BG . (U t H BG U = E where E is the diag- 
onal matrix with elements E n .) The equilibrium density 
matrix is p° = exp(— (3TL)/Z with Z = Tr cxp(—/3H), 
ft = 1/ksT, and the order parameter equation is 



A?, = %[Up°Ut] 



ij+N , 



(2) 



Note that in our notation, half of the quasiparticles states 
are occupied in the ground state. Eq. || is the vortex- 



lattice state analog of the simple Meissner state gap equa- 
tion and can be solved iteratively. 

The linear response of the density matrix at frequency 
u) is given by the following familiar expression: 



f(E n )-f(E n >) 
Tilo — E n i + E n + irj 



[Ut,5H BG LV,„ (3) 



5H BG has a contribution (5e) in its diagonal blocks, 
which describes the coupling of a weak a.c. electric field to 
the system |^l],^2j , and a contribution in the off-diagonal 
blocks, 



SAUu) = V ij \U6p(u;)l^] 



ij+N , 



(4) 



which describes the self-consistent response of the order 
parameter. In non-self-consistent calculations the latter 
contribution is neg lected. Multiplying, Eq. (|) by the 
energy denominantor, leads to a set of linear equations 
for 5p n > n which can be solved by inverting a symplectic 
pl[ matrix that is singular when hoj equals a collective 
excitation energy of the system. Since 



SO = ^2(n'\0\n)8p n,n i 



(5) 



we can calculate the linear response of any observable to 
the perturbing a.c. field. The complex optical conductiv- 
ity tensor a^ v {uj) is defined by SJ^ = af iV (cj)E v (u}), 
where J is the current operator and E(oj) is the a.c. 
electric field. Here we report results for an on-site at- 
tractive pairing model (Vij = — VoSij) which leads to 
s-wave superconductivity. It is convenient to compare 
normal and superconducting state response by adjusting 
the strength of the pairing interaction, rather than field 
strength or temperature. The chemical potential is cho- 
sen near the bottom of the band where the dispersion is 
nearly parabolic. 

Resonances in the conductivity are determined by 
the pole structure of the nonequilibrium density matrix 
5p(oj), shown in Fig. [l[ While 8p(ui) has a complicated 
spectrum, only a small subset of all possible resonances 
are excited by a long wavelength EM field, and the re- 
maining resonances are optically silent. In fact, in the ab- 
sence of pinning and disorder there is essentially only one 
optically active resonance, which occurs at a frequency 
<x>o near u> c , and well below w qp , as shown in the top panel 
of Fig. H where the real part of the circularly polarized 
conductivity [u + (uj) = a xx (u>) + ia X y(uj)\ is plotted. We 
stress that the absence of a peak in the conductivity at 
Wqp is not a product of changes in optical selection rules 
brought about by vortex motion, but rather reflects the 
absence of a pole in 8p(uj) at that frequency, as shown 
explicitely in Figs. [I] (a) and (b). This underscores the 
collective nature of the excitation spectrum of the vortex 
lattice. 

In the absence of extrinsic pinning, the small frequency 
shift luq — lu c of the response is a result of broken trans- 
lational invariance by the lattice. In general, we expect 
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FIG. 1. Spectrum of the nonequilibrium density matrix 
5p(ui) at T = 0. (a) Non-self-consistent calculation, (b) 
Self-consistent calculation, (c) With single pinning center per 
unit cell, (d) With two pinning centers per unit cell. Optically 
active resonances with a large spectral weight are indicated 
by thick solid lines. The weaker spectral weights are smaller 
by at least one order of magnitude. The resonant frequencies 
in (a) correspond to quasiparticle energy differences. Results 
are for a magnetic unit cell containing two vortices and 32 
lattice sites. 



vortex motion to be strongly influenced by the atomic 
lattice when the coherence length is comparable to the 
atomic lattice spacing, a situation which may in fact be 
realised in BSCCO p3[ . The effects of including homo- 
geneous extrinsic pinning are also shown in the top panel 
of Fig. |^. Pinning is modelled here by reducing the pair- 
ing interaction Vo to a smaller value V ' at each vortex 
core center, a change which has no effect on the equilib- 
rium vortex lattice state. Pinning shifts the resonance to 
higher frequencies and reduces the spectral weight of the 
mode, with the missing spectral weight appearing in the 
zero-frequency (superfluid) response (bottom right panel 
of Fig. |^). When V ' — 0, the pinning is strongest and 
the resonant frequency is close to w qp . This is consis- 
tent with expectations: as the order parameter response 
is supressed the collective mode should increasingly re- 
semble quasiparticle pair creation. From the Hall-angle 
sum rule 24,2]]], there is a direct relation between u>q 
and the spectral weight of the mode or, equivalently, the 
penetration depth: 
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where the f-sum rule value, ftot, is proportional to the 
oscillator strength of the collective mode. In Fig. || we 
show the effects of pinning on the zero-temperature pen- 
etration depth. 
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FIG. 2. Real part of the electron-active circularly polarized 
conductivity of vortex lattice states. Top: Vq has the same 
value at all vortex cores. Results of a non-self-consistent cal- 
culation are also shown. Bottom left: Every second vortex 
core is pinned. Bottom right: Effect of pinning on the pen- 
etration depth. Open bars are for inhomogeneous pinning. 
Conductivity is measured in units of tf to t where t is the near- 
est neighbour hopping energy. 



The dependence of penetration depth on pinning 
demonstrates the importance of distinquishing the 
macroscopic penetration depth of a vortex lattice, 
which characterizes its low frequency response to long- 
wavelength electromagnetic radiation, from the local ef- 
fective penetration depth [ p5|]26| ] measured by /^SR or 
NMR experiments p^J . The latter is actually an equilib- 
rium property of the vortex lattice, and an index of the 
inhomogeneity of the internal magnetic field distribution. 
(In London theory, the field distribution is determined 
by the zero-field penetration depth.) The /iSR penetra- 
tion depth is only indirectly related to the low-frequency 
limit of the vortex-lattice conductivity. The distinction is 
most stark when disorder and both extrinsic and intrinsic 
pinning pinning are absent; the macroscopic penetration 
depth should then diverge, while the local /iSR penetra- 
tion depth will stay close to its zero-field value. 

Finally, we address the case of inhomogeneous pinning, 
which we expect to be be typical for experimental sys- 
tems and to differ importantly from the homogeneous 
pinning case typically studied in theoretical models |^,^J . 
Since our magnetic unit cells contain two vortices, we 
have the option of pinning them unequally. Interestingly, 
we find that the conductivity spectrum for the case of one 
pinned vortex per unit cell (bottom left panel, Fig. ^) is 
not congruent with the naive anticipation of independent 
contributions from pinned and free vortices. Instead, the 
collective resonance of the unpinned and homogeneously 
pinned cases is split, with one peak lying at to < uj c and 
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the other at uj c < w < uj qp . If the response were lo- 
cal, a peak would occur near u> c (perhaps blue-shifted by 
inter-vortex forces) and a second smaller peak due to the 
pinned vortex would appear near cj qp . The breakdown 
of the naive picture reflects the highly nonlocal nature of 
the clean vortex lattice response. 

In summary we report on a practical numerical ap- 
proach for the study of equilibrium and linear response 
properties in the mixed state. This work is motivated by 
the growing realization that high T c superconductors are 
in the clean, long quasiparticle mean-free-path limit at 
low temperatures and by the need to avoid simplifying 
approximations which become dubious in this limit. In 
this paper we have examined the influence of pinning ar- 
rangements on the ac conductivity of dense vortex states, 
demonstrating the dominantly collective nature of the re- 
sponse and the essential role of non-locality. This type 
of calculation can be used to test potential paradigms 
and to help develop sound phenomenologies for linear re- 
sponse properties of these complex states. 
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